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ABSTRACT 

We calculate the coefficient of bulk viscosity by considering the non-leptonic weak 
interactions in the cores of hybrid stars with both hyperons and quarks. We first 
determine the dependence of the production rate of neutrons on the reaction rate of 
quarks in the non-leptonic processes, that is r„ = K S T S + Fa + 2IV . The conversion 
rate, K s in our scenario is a complicated function of baryon number density. We also 
consider medium effect of quark matter on bulk viscosity. Using these results, we 
estimate the limiting rotation of the hybrid stars, which may suppress the r-mode 
instability more effectively. Hybrid stars should be the candidates for the extremely 
rapid rotators . 

Key words: dense matter — gravitation — stars: neutron — stars: rotation — stars: 
oscillations 



1 INTRODUCTION 

It has been recognized that weak interaction processes in neutron star matter contribute to the bulk 
viscosity (Sawyer & Soui (11979l) : Sawyer (11989al) : Haensel & Schaeffer (119921) ). When only nucle- 
ons and leptons exist in the interior of neutron stars, the bulk viscosity results from Urea processes. 
Some exotic particles, such as hyperons and quarks, can be produced as the density inside neutron 
stars increases (Weber (2005); Glendenning et al. {1992)). The non-leptonic processes involving 
both hyperons and quarks lead to larger coefficient of bulk viscosity of dense nuclear matter(Jones 



* E-mail:zhxp@phy.ccnu. edu.cn 



2 Na-Na Pan , Xiao-Ping Zheng and Jia-Rong Li 

( |2001a,b| ); Wang & Lu (1984)). Therefore, the calculations of the bulk viscosity induced by non- 
leptonic reactions have received considerable attention in the past few years(Sawyer (1989b); 
Madsen (TT992T) : Goyal et al. (1T9941) : Haensel & Levenfish (12000b : Haensel et al. ([2002]);Lindblom 
& Owen (120021) : Zheng et al. (120031) : Zheng et al. (l2004l : Zheng et al.(|2005]))since the discovery 
by Andersson, Friedman and Morsink of r-mode instability in neutron stars(Andersson (1998); 
Friedman & Morsink (1998 )). Actually, neutron stars can be constituted by a larger variety of par- 
ticles other than just neutrons and protons or pure quark matter. The mixed phases, with nucleon- 
hyperon, nucleon-quark or nucleon-hyperon-quark matter, may be contained in the interior of 
neutron stars(Glendenning (1997); Glendenning (119921) ). The bulk viscosities of mixed nucleon- 
hyperon matter and nucleon-quark matter have been discussed by some authors so far(Lindblom 
& Owen (|2002|) : Drago et al. (|2005|) : Nayyar & Owen (1200510 . However, the calculations are re- 
duced to single reaction case since only one type of reaction is involved, either hadron or quark 
non-leptonic process. The purpose of this paper is to investigate the general case thoroughly. We 
focus on mixed nucleon-hyperon-quark matter to evaluate its bulk viscosity by considering the 
reactions involving both hyperons and quarks. According to Lindblom & Owen (LO) formalism 
(Lindblom & Owen (120021) : Nayyar & Owen (12005|) ). the bulk viscosity depends on the relaxation 
time governed by the hyperon and quark reactions. Following the LO approach, we first convert 
reaction rates for hyperons and quarks into the production rate of neutrons per unit volume in order 
to give the relaxation time of the system, and then simulate the bulk viscosity of the mixed phase 
using the rates. In our study, the conversion rate from quarks to neutrons varies with a variable 
baryon number density, and we will see below that the numerical solution is necessary. Our study 
is under the neutron stars models based on Gho sh- Phatak- S ahu (GPS) equations of state(Ghosh 
et al. (fl995T) ) (EOS) for hadron matter and SGT model (Schertler, Greiner & Thoma (fT997T) ) for 
quark matter in medium. At the high density, GPS EOSs are favorable for the equilibrium state 
to include a sequence of hyperons. Our model also concerns a Baym - Pethick - Sutherland 
(BPS )(Baym et al. (119711) ) crust. Our analysis shows that there exist two instability windows for a 
neutron star with strangeness. The low-temperature window indicates the possibility of existence 
of submillisecond pulsars. 

The organization of the rest of this paper is as follows. In section 2, we provide details of 
structures in stars which we are studying by taking into account a first-order deconfinement phase 
transition into quarks from hadrons in the cores of the stars. Then in section 3, we employ LO 
approach to derive the bulk viscosity of mixed nucleon-hyperon-quark matter, including numerical 
aspects of the evaluations of the viscosity coefficient. Finally in section 4, we work out r-mode 
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instability window based on the competition between increasing of gravitational radiation and 
damping of viscosities. 



2 EQUATION OF STATE 

In this paper, we consider the compact stars inside which the deconfinement transition occurs 
at high density. Concerning the hadronic phase (HP), we use the BPS EOS at subnuclear den- 
sities for the crust of the star, and the GPS EOS for supranuclear density in the framework of 
the relativistic self-consistent mean field theory (Glendenning (1997)). In this theory, E~ and A 
hyperons are included for nuclear densities, where nucleons interact through the nuclear force me- 
diated by the exchange of isoscalar and isovector mesons (cr, co,p). The model parameters used in 
this paper are arranged in table 1. These two EOSs are matched at e « 10 14 g/cm 3 « 0.4eo, here 
6o = l40Mev/fm 3 w 2.5 x 10 14 g/cm 3 is the saturation density of nucleon, above which we allow 
the HP to undergo a first order phase transition to a deconfined quark matter phase (QP) that is 
described with the effective mass bag model considering medium effect represent by the coupling 
constant g for strong action (Schertler, Greiner & Thoma (119971) ). This phase transition makes it 
possible that the occurrence of a mixed hadron-quark phase (MP) in a finite density range inside 
compact stars. The relative constraints of various species in the MP are determined at each density 
by imposing /^-equilibrium, Gibbs condition and global charge neutrality. At the high density of 
interest to us, these equilibrium constraints are 



Hn = Mu + 2//rf, 




(1) 


Hp = 2Hu + fia, 




(2) 


Mfi — Me, 




(3) 


fiX- = Mn + fie, 




(4) 


Ma — Mn, 




(5) 


Ms = Md, 




(6) 


Mu = Md - Me, 




(7) 


PHp(Mn,Me) = P()p(Mn,Me), 


< 2 1 1 

= X(-^ n u + ^ n d + 3^ + n QP,e)- 


(8) 


(1 -x)(n p - n x - - n H p e - n M ) 


(9) 



Here //,-, n { are the chemical potential and the number density of the i species respectively, 
Pi is the pressure of the i tb phase, and x is me fraction of the quark matter in MP. Solving these 
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constraints, we can give the compositions of the equilibrium matter in hybrid stars. Figure [Hand 
[2] show some solutions for different parameters. Evidently, they contain a significant number of 
hyperons and strange quark matter in the highest density portion of their cores, and the additional 
hyperons appear as the coupling among quarks increases. Using the Tolman-Oppenheimer-Volkoff 
(TOV) equation, we obtain the structures of some hybrid stars with maximum masses for a few 
EOSs shown in Figure [31 where we also plot the structures of hyperon stars for the same GPS 
EOSs (Figure (3^c)) as comparison. The EOSs of the stars are softened when deconfinement phase 
transition happens, but the structures of the stars are not sensitive to the changes of GPS EOSs 
unlike hyperon stars, while medium effect in quark matter does significantly change the stellar 
masses and radii. As is known, the Hulse-Taylor pulsar has the measured mass of 1 .44 M mn (Taylor 
& Weisberg (11989I) ). Strictly speaking, the EOSs models producing the maximum mass smaller 
than the measured value should be rejected. We, however, take the lower values around 1 .44 M sun 
to compare with GPS models (hyperon stars), since the small difference of the masses can't change 
our results significantly below. 



3 BULK VISCOSITY 

Bulk viscosity is the dissipative process in which the macroscopic compression or expansion of a 
fluid element is converted to heat(Landau & Lifshitz (1999)). The dissipation due to bulk viscos- 
ity is carried out via microscopic reactions which come from weak interactions in compact star 
matter. In general, there are two types of reactions in dense nuclear matter. One produces bulk 
viscosity because its timescales are comparable with the period of the perturbation, the other puts 
constraints on the variations of different particles' densities. As in most cases, the volume per unit 
mass variation and the chemical potential imbalance due to the oscillation in the star are so small 
that it is the linear part that contributes most to the viscosity, so we will use the relaxation time 
approximation method(Lindblom & Owen (|2002l) : Drago et al. (|20051) : Landau & Lifshitz (| 19991) ) 
to compute the bulk viscosity of dense nuclear matter. 

Here we are interested in the case that hyperons and quarks are all present in the mixed phase. 
It is thought that the non-leptonic processes are dominating over the bulk viscosity of nuclear 
matter with strangeness (Jones fl2001a,bD ; Lindblom & Owen (2002); Drago et al.(2005)). In the 



initial works(Drago et al.(2005 )), one estimated the results arising from a single reaction, assuming 
either an absence of hyperon or quark pairing in mixed hadron-quark phase. That treatment was 
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for simplification in the calculations. In fact, the processes that produce bulk viscosity of the mixed 
system are 

u + d < — > u + s (10) 
for strange quark matter, and 

n + n < — > p + IT, (11) 
n + p < — > p + A (12) 
for hyperon matter. 

Accordingly, the processes which put constraints on the variations of the densities of different 
particles come from the melting reactions of nucleons 

p < — > 2u + d, (13) 
n < — > u + 2d, (14) 
n + A < — > p + I~, (15) 
and hyperons 

A < — > u + d + s (16) 
into free quarks. 

In the mixed phase, since the degrees of freedom, for example, the number densities of various 
baryons, are related to each other even out of thermodynamic equilibrium by constraints such as 
conservation of baryon number, we could express all of the perturbed quantities in terms of a single 
one. We here choose the number density of neutrons n n as our primary variable. Since the three 
reactions (fT0l)([TTI)([T2l) all contribute to the change of neutron numbers, we express the production 
rate of neutrons per unit volume as 

r n = K s Y s + K A T A + K^ , (17) 

where T s , T A , F x - indicate the reaction rates of processes (flOl) . (fTTI) and (fT2l) respectively, and K s , 
Ka, Kx- are the rates converting reaction rates of each exotic particles, s quark, A and £~ hyperon, 
into the production rate of neutrons. For fast processes, we know that their contributions to the 
relaxation time of the system could be omitted compared with the slow processes, so we may 
easily know that K A = I and K x = 2 by the reactions CCD) and (fT2)) . The determination of K s is a 
slightly complicated matter. Its treatment will be given later. Obviously, if the system is perturbed, 
the fraction of the quark matter in the mixed phase x would be changed through reactions (fT3l . 
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(fl4l) and (fT6l) . and the pressure of the system would be equilibrated again after the perturbation in 
order to satisfy the mechanical equilibrium. Therefore we have following constraints: 



= (1 -X)(8n n + 6n p + 8n Y - + 6n A ) + x(8n u + 6n d + Sn s )/3 



+8x((n u + n d + n s )/3 -n n -n p - n x - - n A ), 



(18) 



= (1 -X)(8n p - 6n Y ) + xQ8n u - 8n d - 8n s )/3 




(19) 



(20) 



(21) 



= a nn 6n n + a np 8n p + a nI ,-8nz- + a nA 6n A - a uu 8n u - 2a dd 8n d , 



(22) 



= /3 n 8n n + /3 p 8n p + fc-8nx- + f3 A 8n A , 



(23) 



= a An 8n n + a Ap 8n p + a AI ,-8n Y - + a AA 8n A - a uu 8n u - a dd 8n d - a ss 8n s . 



(24) 



The first constraint ([TBI is the baryon number conservation of the system. The second one (fT9l) 
imposes the conservation of electric charge assuming all leptonic reaction rates are much smaller 
and neglected. The third constraint (|20l) is related to the mechanical equilibrium. The last four 
(|2TI)-(|241) describe the equilibrium with respect to reactions (fT3l) - (fT6l) . Noteworthily, when the 
density vanishes in the system, 8n x - = 0, and the melting process (fT5l) doesn't exist, nor does the 
process (fTT|) . 
Here we use 

a u = {dHildnj) nkMj , (25) 
Pi = a ni + a Ai - a pi - a^-t, (26) 



for shortening. Then we could obtain the variations of various particles' number densities 8n t and 
the quark fraction 8x in the form of 8n n . 

In order to determine K s , we need to describe the variation of neutrons 8n n following ([T7T) as 



Pi = dP/dni 



(27) 



- 8n n = K s 8n s + 8n A + 28n^-. 



(28) 



Thus we have 
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which is a function of baryon number density for given equation of state under thermodynamic 
equilibrium in the mixed phase system. Combining the constraints conditions to equation (l29l) 
with the linear dependence of 6n A , 6n^- and 6n s on 6n n , the K s can be determined. An example 
of the numerical result is showed in the Figure |4] In the extreme, if the system doesn't undergo 
the deconfinement phase transition (6n s = 0) or is in the absence of hyperons (6n A = 6n^- = 0), 
the calculation is reduced to the situation of a single reaction(Lindblom & Owen (120021) : Drago et 
al. (l2005T) ). 

The relaxation time r of the system due to these reactions must be: 

I. + (30) 

t 8ji 8p 5fi 8n n 

Here Y s , Y A and r s - have already been calculated by Dai under the high temperature limit (2nkT » 
8n) (Dai & Lu (fl996l) ) and Lindblom et al. (Lindblom & Owen (120021) : Nayyar & Owen (120031) 1 
separately. Considering the reactions (|T3T)dT4])(fT3T)(fT6b . the overall chemical potential imbalance 
8jj. is obtained through the equilibrium of the system : 

6/x = 6/u d - 6fi s = 6/u„ - 6p A = 2Sp n - 6/u p - Sjj.?-, (31) 

which could also be described in the form of 6n n . 
Now we use the LO formula 

(32 ) 

1 + (.COT) z 

to evaluate the bulk viscosity. Thereinto, 
_ n\ dp dx n 

7co-7o = a , ■ (33) 

p on n dn B 

Here x„ = n n /n B is the relative population of neutrons in the equilibrium state of the mixed phase 
calculated in the part 2. Co is the perturbation frequency in compact stars, which is typically 10 3 ~ 
lO 4 ^ 1 . We have chosen Co = 10 4 5 _1 for calculation . By the equation (|32|) . we know that for 
a perturbation of fixed frequency , the bulk viscosity f is in proportion to the reciprocal of the 
relaxation time r under the low temperature limit 1 <sc Cot , that is f oc r _1 ; on the other hand £ is 
in proportion to r under condition of the high temperature limit 1 » Cor, that is £ oc r. Thereby, 
C,(n B ) could be divided into two different behaviors, in another word, £,(n B ) will increase with the 
temperature for low temperature regime; and the circs is just opposite for the high temperature 
regime. The transition temperature is the function of baryon density and perturbation frequency 
(FigureHfor GPS2,g = 2.0 EOS). 

Figure [5] and Figure |6] display the bulk viscosity of mixed phase as a function of baryon den- 
sity for various temperatures and given EOSs. The numerical results show that i) the bulk viscosity 
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of mixed phase, either nucleon-hyperon phase or nucleon-quark phase or nucleon-hyperon-quark 
phase, is of almost the same order of the bulk viscosity of pure strange quark matter, ii) the dif- 
ferences between EOSs only change the compositions and sizes of the mixed phase in the interior 
of the stars, iii) the transition temperature between low and high temperature regime is about 
10 8 ~ 10 9 K in our models. 



4 R-MODE INSTABILITY WINDOWS OF HYBRID STARS 

The r-modes in a perfect fluid star succumbing to gravitational radiation (GR) drive the Chandrasekhar- 
Friedman-Schutz instability for all rates of stellar rotation, and arise due to the action of Coriolis 
force with positive feedback of the increasing of GR(Andersson & Kokkotas (120011) ; Stergioulas 
(120031) 1. Actually, the viscosity of stellar matter can hold back the growth of the modes effec- 
tively. Based on the competition between the destabilizing effect of GR and the damping effect 
of viscosity, a r-mode instability window can be defined. Now, we could apply our results of 
the bulk viscosity in the previous section together with those of the HP and QP to considering 
the r-mode instability windows of rotating compact stars using the standard formulae given in 
literatures(Lindblom et al.(1998); Lindblom et al. dl999l )). We first integrate the viscosity on the 
structure of stars, which has been showed in Figure [3]for a wide range of EOSs, to obtain the bulk 
viscosity timescale tb 
1 1 dE 

- = ~^(-J7)b, (34) 
t b IE dt 

where 
~ 1 



- -a z Q 2 R 2 e(r)r 6 dr, (35) 
2 Jo 

(^)b = -4tt f £(e(r)) < \3 • c^| 2 > r 2 dr. (36) 
dt Jo 

Here E is the mode energy in a frame rotating with the star, (^)b is the rate at which energy is 
drained from the mode by viscosity, a is a dimensionless amplitude coefficient, Q is the angular 
frequency of a spinning star, e(r) is the energy density of star at radius r, and 5l3 is the Eulerian 
velocity perturbation. In general the expansion of the mode is a complicated function 

of radius and angle. Although the function can only been determined numerically, we adopt an 
excellent analytical fit to those numerical solutions 

< ' >= l^&l 1 + 0-8<) 2 ]& 2 , (37) 
690 R R nGe 
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which is given by Lindblom (Lindblom & Owen (2002); Nayyar & Owen (12005|) ). Finally, we 
could obtain the critical rotation frequency for the star with a solid crust as a function of tempera- 
ture following from 

- — + — + — = (38) 

T GR r B T SR 

for low temperature limit(r < 10 9 ^f), and 

- — + — + — = (39) 

TGR t B TmU 

for high temperature limit(r > 10 9 K). Here r GR is the timescale for gravitational radiation to effect 
the r-mode( Anders son & Kokkotas (2001|) ). while tsr is the timescale for surface rubbing due to 
the presence of the viscous boundary layer if a solid crust exists(Bildsten & Ushomirsky (2000); 
Andersson et al. (|2000l) ). and T m y is the viscosity timescale of modified URCA reactions for low 
density hadronic matter(Lindblom et al.( 1999)). Our results for the critical rotation frequency are 
shown in the Figure[8l 

We may know from it that there are two instability windows for such hybrid stars( Anders son 
& Kokkotas (2001)). We find that the r-mode instability is completely suppressed for the hybrid 
stars from the low temperature windows. We realize that hadronic EOSs, either stiffer of softer, 
hardly influence the instability windows. It is because that the structure of the star including quark 
phase hardly changes for different GPS EOSs. The increasing medium effect of quarks enlarges 
the window at low g values but reduces the window at high g values. The properties may imply 
that extremely rapid rotating or submillisecond pulsars would exist, and just be the hybrid stars. 



5 CONCLUSION 

Hyperons and quarks exist in cores of neutron stars at high density. To evaluate accurately the bulk 
viscosity (due to strangeness-changing four-baryon weak interactions involving both hyperons 
and quarks non-leptonic processes), it is necessary to compute detailed and accurate models for 
the compositions and structures of the neutron-star cores. We use GPS EOSs for hadrons and 
S GT model with medium effect for quarks to evaluate composition of the nuclear matter inside 
the stellar core and solve the TOV equation to obtain the stellar structure, and particularly study 
the effect of differences between EOSs on the structures of stars. 

We numerically calculate the bulk viscosity of mixed phase matter of various EOSs, and es- 
pecially give the treatment of the case that involves both hyperons and quarks reactions. It is a 
generalization of single reaction case in previous work (Drago et al. (120051) ). Our results show that 
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there is high coefficient of bulk viscosity as long as exotic particle appears in high density nuclear 
matter, either strange quark or hyperon. 

We have considered medium effect of quark matter. The bulk viscosity of hybrid star increases 
as the coupling among quarks becomes stronger. The coefficient of bulk viscosity at large coupling 
constant is higher than assumed in reference (Drago et al. (|2Q05| ), and the mixed phase regime is 
expanded. Just so, the hybrid stars have extremely high limiting rotation frequencies from the low 
temperature windows even if the stellar masses are large. The submillisecond pulsars may exist. 

We must emphasize that we are concerned about the bulk viscosity of mixed normal nucleon- 
hyperon-quark phase. Nucleon, hyperon superfluidity and quark superconductivity may occur at 
T < 10 8 A^ and given critical density. These effects would significantly influence the dissipation 
of hybrid star matter. These considerations, which are made in hyperon star model (Lindblom & 
Owen (2002)) and hybrid star model (Drago et al. (|2005|) ). are our future work. 
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Name 


poifm 3 ) 


B/A(Mev) 


a sym (Mev) 


K(Mev) 


m*/m N x a = x p = x 0} 


GPS I 


0.150 


-16.0 


32.5 


250 


0.83 -Jl 


GPS 2 


0.150 


-16.0 


32.5 


300 


0.83 Jf 


GPS 3 


0.150 


-16.0 


32.5 


350 


0.83 J\ 



Table 1. The nuclear matter properties of the GPS equations of state, po is the saturation density, B/A is the binding energy, and the incompress- 
ibility, effective mass and the symmetry energy are denoted by K, m" jm N and a sym respectively. Here we choose hyperon coupling parameters x a , 

x p , x 0) to be times the nucleon-meson coupling parameters. 
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Figure 1. Influences of medium effect on the relative particles populations as functions of baryon number densities in the hybrid stars. Note that 
under the condition g = 1.0, A does also exist in the MP, but its population is small compared with others', so we don't show it here. Here they are 
the relative populations of baryon number densities for u,d,s quark only. For all QP EOSs, we choose the bag constant B l/4 = 170.0Mev and the 
mass of the s quark m s = 150.0M<?v in calculations. 
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GPS1 g=0.0 



GPS2 g=0.0 



GPS3 g=0.0 



0.0 0.2 0.4 0.6 0.8 1 .0 1 .2 1 .4 1 .6 1.8 2.0 



n B (fm ) 



Figure 2. Influences of hadron EOSs on relative particles populations as functions of baryon number densities in the hybrid stars. Under the 
condition GPS 1, A also exists in the MP, but we don't show it here as its population is small compared with others'. Parameters of QP EOSs are 
the same with the Figure 1 in calculations. 
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Figure 3. Structures of maximum masses hybrid stars (a) (b) using different EOSs discussed in the Figure 1 and 2, and the corresponding hyperon 
stars (c) for same hadron EOSs with case (b) : the relations between total energy densities(in units of 10 15 g/cm 3 ) and the distances from the center 
of the stars . 
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Figure 4. Converting rate K s as a function of baryon number density in the mixed phase under the conditions GPS2 and g = 3.0. 
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Figure 5. Influences of medium effect on bulk viscosities as functions of baryon number densities for a range of temperatures. Parameters are the 
same with Figure 1. 
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Figure 6. Influences of hadron EOSs on bulk viscosities as functions of baryon number densities for a range of temperatures. Parameters are the 
same with Figure 2. 



Bulk viscosity of Mixed nucleon-hyperon-quark Matter in Neutron stars 17 



c 

re ies 



1000000 I L 



~i 1 r 



LGPS2g=2^ 




=10 4 s" 1 



&=10 



3 -1 



_i I i_ 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

n B (fm" 3 ) 



Figure 7. Critical transition temperatures as functions of baryon number densities and perturbation frequencies under the conditions GPS 2 and 
g = 2.0. 
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Figure 8. The influences of medium effect and hadron EOSs on r-mode instability windows. Each pair of curves represent the hybrid stars of fixed 
maximum mass discussed in the Figure[5Ia) (b) respectively. Here v is the critical rotation frequency of the stars. 



